Correlated fermions in a one-dimensional quasiperiodic potential 
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We study analytically one-dimensional interacting spinless fermions in a Fibonacci potential. We 
show that the effects of the quasiperiodic modulation are intermediate between those of a com- 
mensurate potential and a disordered one. The system exhibits a metal-insulator transition whose 
position depends both on the strength of the correlations and on the position of the Fermi level. 
Consequently, the conductivity displays a power law like size and frequency behaviour characterized 
by a non trivial exponent. 
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Since the discovery of quasicrystals (QCs) in 1984 
by Schechtman et al. |5J], the electronic properties of 
quasiperiodic systems have been intensively studied. 
These metallic alloys are notably characterized by a 
low electrical conductivity a which increases when ei- 
ther temperature or disorder increases H. The very 
low temperature behaviour of a is still an open ques- 
tion and depends on the materials. For example, in Al- 
CuFe H and AlCuRu a finite conductivity at zero 
temperature is expected whereas recent results seem 
to confirm a Mott's variable range hopping mechanism 
0(T) - exp -(To/T) 1 / 4 ) for i-AlPdRe icosahedral phase 
down to 20 mK. The optical conductivity is also unusual 
since there is no Drude Peak for icosahedral quasicrystals 



Many theoretical works have attempted to understand 
how the quasiperiodic order could induce such exotic be- 
haviours. In particular, the case of independent elec- 
trons in one-dimensional (ID) systems has been deeply 
investigated for different structures (Harper model, Fi- 
bonacci chain,..) |^,^|, giving rise to singular continu- 
ous spectra with an infinite number of gaps. Moreover, 
the corresponding eigenstates are neither extended nor 
localized but critical, and are known to be responsi- 
ble of anomalous diffusion ]lO| , |ll| . For higher dimen- 
sional systems (Penrose tiling, Octagonal tiling, Icosahe- 
dral structure...), similar studies had also displayed com- 
plex and intricated spectra, with analogous characteris- 
tics of the electronic states [[l2],[l3|. The transport prop- 
erties have also been widely studied |L4]-[L8| within differ- 
ent approach (Landauer conductance, Bloch-Boltzmann 
approximation, Kubo-Greenwood conductivity,...). How- 
ever, given the complexity of these problems due to the 
geometry alone, the interactions between electrons have 
often been neglected. Even in ID incommensurate struc- 
tures, few results have been obtained [p^|-p3[. 

In this Letter, we investigate the effect of the inter- 
actions considering a Hubbard-like model for spinless 



fermions embedded in a Fibonacci potential. We take the 
correlations into account using a bosonization technique 
whereas the quasiperiodicity is treated perturbatively. 
Using a renormalization group approach, we show that 
the quasicrystalline system displays a metal-insulator 
transition (MIT) induced by the interactions. The cor- 
responding critical MIT point is found to be strongly 
dependent on the Fermi level. In marked contrast with 
the simple cases of disordered or commensurate poten- 
tials, the quasiperiodicity leads to a power law like de- 
pendence of the conductivity either in size or frequency 
with an exponent depending both on the interactions and 
on the position of the Fermi level. Though our analysis is 
performed for a Fibonacci potential, we stress that these 
results can be extended to any potential having a non 
flat dense Fourier spectrum. 

Let us consider a model of interacting spinless fermions 
on a quasiperiodic lattice described by the following 
Hamiltonian: 



rii rij 
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(i) 



where c] (resp. c;) denotes the creation (resp. annihila- 
tion) fermion operator, rii = cj Ci represents the fermion 
density on site i, and (...) stands for nearest neighbors 
pairs. The quasiperiodicity is provided by the Wi's that 
take two discrete values Wa = +A/2 or Wb = —A/2 
given by the spatial modulation of the Fibonacci chain. 
In fact, we consider a periodic approximant of this struc- 
ture with Fi sites per unit cell that can be obtained by I 
iterations of the substitution rules: A — > AB, B — > A, 
where Ft is the I th element of the Fibonacci sequence 
defined by: 



Fi 
Fi+i 
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= 1 
Fi-i 



(2) 



We denote p = F/_ 2 , s = Fi-\,n = Fi and n' — s (resp. 



1 



n' = p) if I is even (resp. odd). In the quasiperiodic 
limit (/ — > oo), the ratio s/p converges toward the golden 



mean r 



It is also useful to compute the Fourier 



transform of the potential W. This can be done via the 
conumbering scheme p4|: 
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(3) 



for to = 1 to n — 1 (a is the lattice spacing). A global 
shift of the Wi allows us to deal with a zero-averaged 
potential so that we can set W(0) = 0. The interest 
of this potential lies in the fact that it is perturbatively 
non trivial (the Fourier transform is dense in [0, 2ir], in 
the quasiperiodic limit). Moreover, the underlying sub- 
stitution rule provides a self-similar structure that can 
be readily seen in Fig. |l|. 



FIG. 1. Fourier transform of the diagonal Fibonacci po- 



tential for the 15 
cell. 



approximant with F15 — 610 sites per 



If the quasiperiodic perturbation is small enough, we 
can focus on the low energy properties around the Fermi 
level. To treat the interactions, it is convenient to use a 
bosonic representation of the fermion operators. In this 
representation the hamiltonian (JIJ) becomes: 

H = H Q + H W (4) 

where H stems from the periodic part (A = 0) of ([!]) and 
reads in the continuum limit [E5|: 



expressed in terms of the microscopic parameters t and 
V: 



uK = vf 



K 



vf 
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7T 



(G) 



Actually, the representation (||) is more general [^6 27 
and gives the correct low energy description of the sys- 
tem, even when the interactions are strong provided the 
exact u and K parameters are used. The quasiperiodic 
part Hyy can also be written in terms of the boson fields 
& 



Hw = 



1 

2na 



dx W(x) cos [2k F x + 2(j)(x)} 



(7) 



where kp is the Fermi wave vector and a is a short dis- 
tance cut-off of the order of the lattice constant a. 

The effect of the quasiperiodic potential is computed 
using a perturbative renormalization group (RG) ap- 
proach, similar to the one for a single harmonic p9| . 
The RG equations for the potential and the interaction 
parameter read: 



dK 
~dT 



= -K 2 G(l) 



G ^= E E^ i? t('? + ^fcF)a(/)] 
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(10) 
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where y q — aW(q)/u is the dimensionless Fourier com- 
ponents of W, a(l) = a(0) e l is the renormalized short 
distance cutoff. In (|l^), the sum over q is performed for 
q = 2irm/n with m € [l,n — 1], and R is an ultravio- 
let regulator whose precise form depends on the cut-off 
procedure used. Without loss of generality, we choose a 
gaussian regulator R(x) = e~ x . Note that the renormal- 
ization of the velocity u is neglected since it is of higher 
order in the potential amplitude A. 

The physical properties of the system are determined 
by the long distance behaviour of K. If K converges to- 
ward a fixed point, the system remains metallic but with 
renormalized parameter (K*,u*). Otherwise, the solu- 
tion flows to a strong coupling regime whose physics de- 
pends on the precise nature of the potential W . Note that 
expression (|^) is valid for any lattice potential perturba- 



(w!f)(7rn) 2 + ( — ) (d x (/)) 2 . (5) tively treated. For a single harmonic (W(q) 



In ([3]) (j) is a boson field related to the long wave length 
part of the fermionic density by p(x) = —'V4>(x)/w, and 
LI is its canonically conjugate field. All the interactions 
are absorbed in the two constants u and K, where u is the 
renormalized Fermi velocity (in the non-interacting case 
(K = 1), u = vf = 2tosin(fci?a)), and K is the parame- 
ter controlling the decay of various correlation functions. 
For weak interactions, u and K can be perturbatively 



two situations occurs. If go 7^ 2kp, R stops the renormal- 
ization of K at large enough lengthscale and the potential 
is irrelevant. A contrario, if go = 2fcj? one recovers the 
usual Metal-Insulator transition at K c = 2. In the disor- 
dered case, the W(g)'s are given by an uniform averaged 

distribution: W* (q)W(q / ) = yS qq r. In the limit of weak 
disorder, (||) can be integrated neglecting the renormal- 
ization of K: y(l) = y(Q)e^ K ^ 1 . Then, Eq.(|) simply 
becomes: 
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^ = -K 2 Ce^ 2K ^ L (11) 

where C is a constant. Eg. (pi]) defines a critical value 
K c = 3/2, separating an insulating phase (K < K c ) from 
a metallic state (K > K c ) |3(J. 

We now investigate the case of the Fibonacci potential. 
For a given maximum renormalization length scale l ma x 
corresponding to accessible physical range, different cases 
must be distinguished. First, if the Fermi momentum 
2k p is close to a main peak of the Fourier spectrum (^), 
then, at long distance (i. e. I ~ l m ax), the flow of K is 
controlled by this harmonic (see Fig. ^ (upper curve)) 
and its behaviour is similar to the periodic one. 



asymptotically correct or not, one would need an analyt- 
ical calculation of G, a rather complicated task. In this 
context, the flow of K is given by: 

^£ = -K 2 De( 4 - 2K -^ 1 (13) 
where D is a constant. 



FIG. 3. Behavior of G for 2&f = 7r (lower curve) and 
2kF = 1.95 (upper curve). Both cases have been offset for 
clarity. 



FIG. 2. Behavior of G for 2/cf = 2.63 (lower curve) and 
2fcp = 2.4 (upper curve). 

There is a transition at K c = 2, with a metallic phase 
for K > 2 where the quasiperiodic potential is irrelevant. 
In this phase, the system has gapless excitations whose 
correlations are given by (^J) with renormalized parame- 
ters. For K < 2, the quasiperiodicity is relevant and the 
system has a charge gap given by: 



V2k f 



(12) 



Note that for K = 1 (non-interacting case), one obtains 
a linear scaling of the gap opening. Our methods allows 
us to recover very simply the perturbative results ]3lf| de- 
rived by quite different methods for the particular case of 
non-interacting electrons. The effect of the interactions 
is thus essentially to change the scaling of the gaps and 
to allow a MIT at K c — 2 (attractive interactions). 

A more unusual behavior is encountered when the 
Fermi level is far from a dominant harmonic of the 
quasiperiodic potential. Indeed, the low energy prop- 
erties up to l max are no more dominated by the ultimate 
presence of a gap or not but by the precise dependence 
of G with the scale. The specific feature of the quasiperi- 
odic case is that contrary to the disordered case, G a 
priori depends on the Fermi level. As shown in Fig. |^, 
it is reasonable to approximate this behavior by an ex- 
ponential scaling. To know whether such a description is 



The position of the MIT point is then given by K c = 
2 — /i/2. In addition, the Fibonacci potential seems to 
provide a unique value ju ~ 2 (see Fig. ||), that leads to a 
transition for the non-interacting point K c = 1. It would 
be interesting to know if this result remains true beyond 
the perturbative theory, and if it is a generic property 
of self-similar potentials. Finally, note that intermediate 
cases can also occur (see Fig. || (lower curve)) for which 
G can not be naively approximated by an exponential 
law. In this context, one can not simply extract a critical 
behaviour from the RG equations. This deserves further 
investigations. 

So, the quasicrystal differs from a periodic one, for 
which the gap only acts for a given position of the Fermi 
level with K c = 2, and a disordered system for which 
the potential is relevant regardless of the position of the 
Fermi level, but below a constant critical value K c = 3/2. 
This important modification of K c is reminiscent of a 
correlated disorder with long range correlations in space 

for which the averaged disorder potential W* (q)W(q') = 
8qq'A(q) is not constant. 

The consequences of the scale dependence of G can 
be directly seen on transport quantities. In the regime 
where the RG equations are valid, the conductivity can 
be computed using a perturbative method. The simplest 
is to use the so-called memory function formalism which 



is well suited to the bosonization representation 29 1 . We 
refer the readers to |32| for the technical details and only 
give here the main results. The high frequency optical 
conductivity and the dependence of the resistance with 
respect to the system size L = e l are given by: 



3 



R{L) oc LM(l L ) 



(14) 
(15) 



where in each case the memory function M is computed 
at a renormalized frequency scale l u = log(t/w) or at a 
renormalized size scale It, = log(L/a(0)). M is pertur- 
batively given by |33| : 



M(l) oc G(l)e 



(16) 



This expression remains valid as long as one remains in 
the perturbative regime, i.e. as long as the renormaliza- 
tion of K in (||) remains small. This is true for all length 
scales in the metallic regime K > K c , where the potential 
W is irrelevant. For K < K c this holds until a certain 
length scale Lsc corresponding to the strong coupling 
limit. Using ( JIB] ) and ([l3]) one straightforwardly gets the 
resistance for the quasiperiodic system: 



R(L) ~ L 



4-2K- 



(17) 



so that R(L) ~ L 4 - 2K , R(L) ~ L 3 ~ 2K for the commen- 
surate and disordered cases respectively. As foreseen, 
in the regime where the scattering potential is relevant 
the resistance increases with the scale. For noninteract- 
ing (K — 1) disordered electrons we simply recover the 
Ohm's law, which is valid below the localization length. 
For a quasicrystal with a dense Fourier spectrum, we find 
a non universal power law increase of the resistance, with 
an exponent depending both on the interactions and on 
the position of the Fermi level. Correspondingly, there 
is a power law frequency dependence of the optical con- 
ductivity <j{uj) ~ (1/lu) 5 ~ 2K ~ 1 - 1 . The properties of po- 
tential W have thus a direct impact on the increase of 
the resistance. This is in contrast with both the disor- 
dered and the commensurate cases where the scaling of 
the resistance only depends upon the interactions. Such 
a behavior and the relationship between the conductiv- 
ity exponent and the spectrum of the quasiperiodic po- 
tential should be testable in numerical simulations. For 
the noninteracting case, some numerical observations of 
a power law behavior for the conductivity have already 
been reported |16|. However the system sizes considered 
were too small to allow a precise determination of the 
exponents, so that more numerical studies are obviously 
needed @. 

For K < K c the previous perturbative analysis ceases 
to be valid beyond the length scale Lsc for which 
G(lsc) ~ 1- For the commensurate and the disorder case 
this defines the correlation length respectively associated 
to the gap and to the localization length. Above Lsc the 
resistance grows exponentially with the size of the system 



R(L) 



,L/L S 



For the quasiperiodic case our analysis 



shows that a corresponding typical length exists as well. 
However its physical interpretation and the physical be- 
havior above Lsc are still to be understood. Since for a 



non-interacting system the wavefunctions of a quasiperi- 
odic system exhibit algebraic decay, a reasonable guess 
is that above Lsc the resistance keeps on increasing as 
a power law but very likely with another exponent than 
4 — 2K — /j. Clearly the investigation of this strong cou- 
pling regime deserves further studies. 

We thank R. Mosseri, B. Doucot and CI. Aslangul for 
fruitful discussions. 
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